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L'invention conceme un precede d'^valuation de la 
degradation tfune image video cod6e par blocs de points 
images ou pixels, li comporte las etapes sulvantes: 

a) s^lectionner une image d'entree et determiner son ac- 
tivlte spatlale (SA) dans une fenetre d'analyse selon les 
sous-6tapes suivantes: 

I) detenminer pour chaque bloc (n, m) les coefficients 
transformes Fn, m (i, y) par une transformee par blocs 

ii) detemniner a partir des coefficients transformes Fn. m 
(i, y) t'actlvite spatlale bsa de chaque bloc 

til) d^temniner I'actlvlte spatlale globale SA^ de la fen§tre 
tfanalyse 

b) selectionner I'image de sortie conrespondant a Tima- 
ge d'entree et determiner ladite activite spatiale globale SAo 
de rimage de sortie 

c) comparer les activites spatiales (SAg) et (SA^). 
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PROCEDE D- EVALUATION DE lA DEGRADATION D'UNE IMAGE VIDEO 
INTRODUITE PAR UN SYSTEMS DE CODAGE ET/OU DE STOCKAGB 
ET/OU DE TRANSMISSION NUMERIQUE 

La present e invention a pour objet un precede 
d'evaluation de la degradation d'une image video, 
mtroduxte par un systeme de codage et/ou de stockage 
et/ou de transmission numerique, notamment un systeme 
mettant en oeuvre un codage et/ou un stockage et/ou une 
transmission de signal video numerique avec un has debit. 

Les degradations sont en general dues aux 
erreurs introduites par les algorithmes de reduction de 
debxt et/ou par les liaisons defectueuses de transmission 
ou bien encore aux defauts des codeurs et des d6codeurs. 

La num6risation des signaux video .a ouvert la 
possibilite de pouvoir copier, stocker ou transmettre 
I'xnformation video en maintenant une qualite constante. 

Cependant, la grande quantite d • information 
vehxculee par les images video n^cessite en pratique 
1' utilisation de methodes de compression numerique pour 
redduire le debit binaire. 

Une methode de compression tres repandue dans 
le domaine de la video est decrite dans la norma ISO/CEI 
13918 MPEG2. Get algorithme est dit du rype "avec 
pertes", car I'image restitute apres le decodaae n ' est 
pas identique a 1' original. Afin de maintenir une aualite 
acceptable pour le telespectateur final, les algorithmes 
de reduction de d^bit tiennent compre des oroprietes 
perceptuelles du systeme humain de vision. En depit de 
ceci, les cont raintes imposees de debit ou de largeur de 
bande disponible pour la transmission ainsi aue le 
contenu du signal video impliquent 1' apparition de 
degraaations caracteristiques sur le signal apres 
decodage. Ces degradations introduites par la chaine, par 
exempie MPEG2, influent directement sur la qualite pergue 
de I'image finale. 
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L ' evaluation automatique de la qualite des 
signaux audiovisuels a un large champ d ' applications dans 
la chaine de television numerique : production, 
distribution, et evaluation des performances des 
5 systemes- 

Par ailleurs, les dispositifs existants ont 
ete elabores pour des tests en laboratoires et ne sont 
pas adaptes pour la telesurveillance des reseaux de 
distribution. 

Des sequences d' images codees MPEG et 
transmises a bas debit pour la diffusion de la television 
numerique ou pour d'autres applications multimedia, vont 
presenter un certain nombre de defauts ou de deformations 
par rapport a la sequence originale. En effet, une liste 
15 non exhaustive de degradations visibles peut etre 
dressee. Les plus perceptibles d'entre elles sont les 
erreurs granulaires, les deformations de contours, les 
pertes d ' information, les contours "exotiques", les 
effets de bloc, etc.. Cependant, des erreurs de 
transmission peu importantes peuvent se traduire par des 
effets plus ou moins localises sur 1 ' image . En cas de 
perturbation imporrante, elles peuvent se traduire par 
des difficultes d'acces aux inf orr:ia-cions , par exemple 
coupures de service cu gels d' images pendant un temps 
plus ou moins long selon la perturbation. L'envergure des 
erreurs depend du niveau de pertinence et de la structure 
des donnees qu' elles affectent : mots de synchronisation, 
vecteurs mouvement, images codees avec ou sans prediction 
ou. images de base pour les predictions. En plus des 
coupures ou gels d* images, les degradations constatees se 
materialisent par des blocs ou des macro-blocs errones ou 
mai posit ionnes. Ceci a pour effet de propager les 
degradations sur toute la sequence video jusqu'a 1 ' image 
codee sans prediction, done codee independamment des 
35 autres. 
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Un precede d' evaluation a ete propose par le 
NTIA (National Teieconununications and Information 
Administration) dans 1 'article de A. A. WEBSTER et 
Collaborateurs, intitule "An objective videoquality 
assesment system based on human perception" et publie en 
juin 1993 dans la revue SPIE vol- 13, p. 15-26. 

Ce precede met en oeuvre une analyse des 
images degradees et des images d'origine, apres les avoir 
filtrees par deux operateurs" de SOBEL (matrices 3x3) 
vertical et horizontal. L' image filtree est obtenue par 
convolution en faisant glisser les matrices de SOBEL 
horizontalement et verticalement , et les resultats 
obtenus representent les gradients verticaux et 
horizontaux de 1* image. En d'autres termes, 1 • image 
15 filtree met en valeur les contours verticaux et 
horizontaux contenus dans 1 'image initiale (non filtree). 

Une mesure basee sur cette information permet 
de mettre en evidence le changement de contenu entre 
1' entree du systeme video et sa sortie. 

Le procede propose par le NTIA met en oeuvre 
deux parametres : 

- d'une part, 1 ' information spariale SI qui 
represente I'ecart moyen mesure sur les pixels de i ■ image 
filtree par I'operateur de SOBEL. II s'agit ici de 
25 determiner I'ecart moyen au niveau des contours de 
1' image filtree, en considerant que les contours sont 
importants pour la vision et qu'ils sont affectes par les 
differents traite.menrs des sysremes numeriques a 
reduction de debit ; 

d' autre part, 1 ' information temporelle TI 
qui represente I'ecart-type de 1' image de difference 
entre deux images successives, cet ecart-type etant 
calcule a partir des differences enrre les valeurs des 
memes pixels de deux trames successives. Le parametre TI 
peut reveler un mouvement saccade du ^ un defaut du 
codeur . 
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Le precede propose par le NTIA met en oeuvre 
un calcul compararif des parametres SI en TI sur le 
signal video numerique, enrre une i.T.age d* entree en une 
image de sortie d'un systeme. 
5 Ce procede presence un certain nombre 

d' inconvenient s . 

Le filtre de SOBEL conserve seulement 
certaines frequences de 1 ' image qui permettent de prendre 
en compte la perte de contour; ce qui fait que la perte 
10 de definition qui pourrait etre prise en compte est ceile 
qui se situe dans la gamme des frequences conservees. En 
d'autres termes, la perte de definition ne peut etre 
prise en compte que de maniere partielle. 

D* autre part, le parametre SI prend en compte 
15 des defauts qui tendent a se compenser. En effet, la 
perte d ' information de 1 ' image tend a faire diminuer le 
parametre SI alors que les faux contours et les effets de 
bloc tendent au contraire a le faire augmenter, ce qui 
fait que le parametre SI n'est significatif que si I'un 
20 ou 1' autre des phenomenes est dominant. 

Enfin, la methode de calcul des parametres SI 
et TI, en utilisanr un ecart type calcule de maniere 
globale sur 1 ' image entiere, reduir de maniere drasnique 
1 ' impac - sur ces parametres des degradations localisees. 

La presente invention a pour objet un procede 
qui permet de remedier au moins en parrie aux defaurs 
precit6s , 

Le procede selon la presente invention 
exploire une transformee par blocs, par exemple la 

30 transformee en cosinus discrete utilisee notammenr dans 
le standard MPEG, pour mettre en valeur des signatures 
caraczeristiques des defauts identifies. Cette approche 
originale permet non seulemen:: de mesurer finement les 
erreurs introduites, mais elle permet egalement de tenir 

35 compte du contenu initial du signal video et des 
algorithmes mis en oeuvre dans MPEG. 
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Les transformations par blocs d*une image 
(transformee de Fourier, transformee en cosinus discrete 
TCD, etc..) sont obtenues par 1* operation : 

[Fn,m] = [T] . [fn.ml . [T]"" ou f (x,y) designe le 
bloc d' image a transformer et T(x,y) la matrice de la 
transformation. One autre transformee par blocs est 
realisee a partir de la transformee en ondelettes de 
1' image, en reorganisant les coefficients d' ondelettes de 
maniere a obtenir des blocs transformes ayant la taille 
desiree et en particulier de meme taille que les blocs 
obtenus par les methodes precitees . Un tel mode de 
reorganisation est indique dans 1" article de 
R. de QUEIROZ et Collaborateurs , intitule "Wavelet 
transforms in a JPEG - like Image Coder", publie en avril 
1997 dans la revue IEEE Transactions On Circuits and 
Systems for Video Technology, vol. 7 n^'Z p. 419-424. 

L'idee de base de 1' invention est en 
particulier de r6aliser des calculs sur les blocs selon 
lesquels 1 ' image transmise a ete codee^ de maniere a 
generer un parametre significatif exempt de I'effet de 
bloc. 

L' invention concerne ainsi un precede 
d' evaluation de la degradation d'une image video codee 
par blocs de points images ou pixels, cette degradation 
etant generee par un systeme de codage et/ou de stockage 
et/ou de transmission produisant une image de sortie a 
partir d'une image d' entree, caracterise en ce qu'il 
comporte les etapes suivantes : 

a) selectionner une image d' entree et 
determiner une activite spatiale (SA) de 1 ' image d' entree 
dans une fenetre d' analyse representant au moins une 
partie de 1 • image qui presente un ensemble de dits blocs 
de points images ou pixels, cette determination mettant 
en oeuvre les sous-etapes suivantes : 

i) determiner pour chaque bloc (n,m) du dit 
ensemble de blocs de pixels, les coefficients transformes 



3NSDOaD: <FR_27e945aA1J_> 



2769453 



6 

Fn,m (i,j) par une transformee par blocs selon les dits 
blocs 

ii ) determiner a parrir des coefficients 
transformes Fn,m (i^j), I'activite spatiale bsa de chaque 

5 bloc dudit ensemble de blocs de pixels 

iii ) determiner a partir de I'activite 
spatiale de bsa de chaque bloc l*activite spatiale 
globale SAx de 1 ' ensemble des blocs de pixels constituant 
la fenetre d* analyse 

1^ b) selectionner 1' image de sortie 

correspondant a 1 ' image d' entree et determiner ladite 
activite spatiale globale SA2 de 1 ' image de sortie dans 
ladite fenetre d' analyse, en mettant en oeuvre les sous- 
etapes a i) ^ a iii) 

^5 c) comparer I'activite spatiale globale (SA2) 

de 1 ' image de sortie dans la fenetre d* analyse et 
I'activite spatiale globale (SAj) de 1' image d' entree dans 
la fenetre d' analyse. 

L' activite spatiale bsa d"un bloc (n,m) peut 

20 etre determinee par la formule suivante : 



35 



L' activite spatiale bsa peut etre obtenue 
comme indique ci-dessus par combinaison des carres des 
composantes de 1 ' image a partir de la quasi-totalite des 
composantes de la TCD. Toute autre fonction de 
30 combinaisqn des composantes de la transformee utilisee 
est applicable pour caracteriser le contenu de 1' image. 
Une fonction plus generale est : 

V'^y=0 ) 



ou k(i,j) est un coefficient constant de ponderation de 
la composante i,j utilisee, et p une constante. 
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II y a plusieurs possibilites pour choisir les 
constanres k(i,j), et: le choix de I'un ou I'aurre se fait 
en foncrion de 1 ' application recherchee. En effet, la 
fonction k{i,j) est choisie suivanr la caracteristique a 
5 mettre en exergue : tenir compte du systeme visuel ou 
exrraire une partie de 1 ' information pertinente. Les cas 
suivants sont proposes : 

1) le procede numero un est a utiliser lorsque 
I'on privilegie une bonne correlation avec la perception 

10 humaine. 

Les valeurs des constantes k(i,j) sont 
initialisees en fonction de 1 ' importance relative des 
coefficients de la transformee pour I'oeil humain, 
notamment sa sensibilite f requentielle, afin de fournir 

15 un parametre d'activite representatif de ce qui est 
pergu. Par exemple, dans le cas de la TCD, on prend 
k(i,j) = l/0(i,j). Les Q(i,j) sonr les composantes de la 
matrice de quantification exploitee pour la reduction de 
debit, indiquees dans le document suivant qui est extrait 

20 de la norme MPEG-2 : ISO/IEC CD 13818-2 : "Information 
technology - Generic coding of moving pictures and 
associated audio information - Part 2 : video", Dec 1993, 
p 45, § 6-3.7. 

2) le second procede esr utilise lorsque les 
25- coefficients par exemple TCD affectes par la degradation 

sont identif iables - Les constantes k(i,j) sonc choisies 
de maniere a eliminer cemains coefficients de la 
transformee utilisee. En effet, la ponderation k(i,j) est 
utilisee pour retenir les coefficients les plus affectes 

30 ou sensibles a une degradation donnee. II s'agit d'une 
matrice binaire affectant zero aux coefficients a 
eliminer et un aux coefficients pertinents. La selection 
se base soit (a) sur la position du coefficient dans la 
matrice par exemple TCD, soit (b) sur son amplitude 

35 moyenne : 
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a) les coefficients correspondant aux hautes 
frequences spatiales sont souvent ies plus affecres par 
la compression. Un exemple de matrice de ponderation 
suivant: le rang du coefficient par exemple TCD est donne 
5 dans ie tableau ci dessous : 

00001111 
00001111 
0 0 0*0 1111 
^0 00001111 

11111111 
11111111 
11111111 
11111111 
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b) certains coefficients de faibie amplitude 
sont ramenes a zero lors de I'etape de compression. 

Afin de choisir ces coefficients^ une moyenne 
ponderee de chacun des coefficients sur la zone d' image 
analysee (de taille M.N. blocs) est effectuee : 

mo,Coe/XW)= I il^^T^ .oCaQ{i,j) est defini 
m=in=il Q(jjJ) \ comme ci-dessus, en 1) 

Les coefficients dont les moyennes sont parmi 
les 48 plus faibles valeurs sont retenus. Pour ces 
derniers, ic(i,j)=i et pour les autres k(i,j)=0. 

Dans le cas du parametre "activite spatiale", 
on fixe k(i,j)=l, et p=2 pour la suite du document 
relative a la description des figures. 

L' activite spatiale globale SA de 1' ensemble 
des blocs peut etre alors determinee par la formule 
suivante : nr., 



5>f = 



osa^ 



35 HxW . 

(HxW) representant le nombre de blocs de 
pixels dans la fenetre d' analyse . 
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Ladite comparaison (etape c) est 

avancageusemeni: effectuee a I'aide du parametre LR defini 
de la maniere suivante : 

LR = gj[fi(SAi,SA2)] avec i,j^{1^2} 
avec /,(x,v)=x-v ou />(x,y)=.r/y 
et g,(x)=100.|x| ou ^,(.T)=IOO.|log(|.r|)| 
et par exemple : LR = logio (SAe/SAs) . 

Le procede peur etre alors caract6rise en ce 
qu'il met en oeuvre le calcul du parametre LR sur une 
sequence de M images d' entree et de M images de sorties 
correspondants, et en ce qu'il met en oeuvre le calcul 
d'un parametre d ' evalua-cion ou note de qualite MLR defini 
de la maniere suivante : 

MLR = Maximum^ (LR) 

Le procede peur etre avantageusement 
caracterise en ce qu'il comporte la determination d'un 
indicateur d'activite temporelle TA pour un groupe de M 
images, determine de la maniere suivante : 
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esA(u) , pour u variant de 0 a M-1, designant 
les iM coefficients d'une transformee par blocs, par 
exemple une transformee en cosinus discrete appliquee a M 
echantillons temporels de I'activite spatiale SA. 

Le procede peut etre caract^ris^ en ce qu'il 
compor-e, en vue de la determination de I'effet de bloc, 
les etapes suivantes : 

determiner I'activite spatiale SA^i de 
1' image d' entree dans une fenetre d' analyse decalee 
presentant des blocs de points image ou pixels qui sont 
decales d'au moins un pixel dans la direction des lignes 
de 1* image et/ou dans une direction perpendiculaire a la 
35 direction des lignes de 1' image par rapport auxdits blocs 
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de pixels de ladite fenetre d' analyse, cette 
determination mettant en oeuvre les sous-etapes suivantes 

i) appliquer ladite transformee par blocs sur 
chaque bloc decale (n,m) dudit ensemble de blocs de 

5 pixels de la fenetre d' analyse decaiee pour determiner 
les coeeficients transformes Fon.m (i,j), 

ii) determiner a partir des coefficients 
transformes F^rxria (i,j) des blocs decales, I'activite 
spatiale bsa^j de chaque bloc dudit ensemble de blocs 

10 decales, 

iii) determiner I'activite spatiale globale 
SA^i de 1* ensemble des blocs decales constituant la 
fenetre d' analyse decaiee, 

e) determiner I'activite spatiale globale SAd2 
15 de 1* image de sortie dans ladite fenetre d* analyse 

decaiee en mettant en oeuvre pour 1 ' image de sortie les 
sous-etapes d i) a d iii) 

f) comparer d'une part I'activite spatiale 
globale SA^2 de 1 ' image de sortie dans la fenetre 

20 d' analyse decaiee et I'activite spatiale globale SA2 de 
1' image de sortie dans la fenetre d' analyse, pour evaluer 
I'effet de bloc dans 1 • image de sortie er d'aucre part 
I'activite spatiale SA^i de 1 ' image d' entree dans la 
fenetre d* analyse decaiee et I'activite spatiale SA^ de 

25 1' image d' entree dans la fenetre d' analyse pour evaluer 
I'effet de bloc dans 1' image d' entree. 

L'effet de blocs se caracterise de deux 
manieres differentes : 

intrinsequement, c'est-a-dire par une 

30 fonction representative du contenu de 1 ' image et de 
1' image seule. 

SM2^gj[MSAj2ySA.)] 

- de maniere dif f erenciee, c'est-a-dire qu'il 
35 est necessaire de comparer deux valeurs d'une fonction 
representative du contenu de 1 ' image : la premiere 
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calculee sur 1 ' image de reference (image d'enrree) et la 
seconde sur 1 ' image degradee (image de sortie) . 

3M = g,[/^{f,{SA^.,SA.),/XSA^^,SA,))], sivec iJ,ke{U2} 

5 

et on a /, (j:, v)= .x- v ou /.{x,y)=x/y 

^,(.r)=100.|.t| ou ^2(x)=I00.|log(|xl)| 

En particulier, la formuie BMj = 100.'^^^^^ est 

10 utilisee. 

L'etape f met avantageusement en oeuvre la 
determination d'un indicateur BM d'effet de bloc, selon 
la formuie suivante : 

SAd2 / SA2 

15 BM = X 100 

SAdi / SAi 

D'autres caracteristiques et avantages de 
1' invention apparaitront mieux a la lecture de la 
description qui va suivre, donnee a titre d'exemple non 
Ixmitatif, en liaison avec les dessins ci-annexes, dans 
lesquels : 

la figure 1 est une representation 
sch6matique d'un dispositif pour ia mise en oeuvre du 
procede selon la presente invention, 

la figure 2 represente une grille alignee 
pour le calcul du parametre d' Evaluation SA et une grille 
decalee pour le calcul du parametre 

la figure 3 represente 1' evolution du 
parametre BM au cours d'une sequence film6e particuliere, 

- et les figures 4a a 4c representent 
respectivement 1' evolution temporelie du parametre SA 
(4a) et la TCD appliquee a M echantillons successifs de 
SA (4b et 4c) . 

Selon la figure 1, des i.-riages d' entree 1^ sont 
codees par un codeur COD et transmises par un systeme de 
transmission TRANS, par exempie un systeme de 



20 
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transmission par satellite, puis decodees a la reception 
par un decodeur DSC pour donner une image de sortie 
decodee I3. Un circuit TCDE applique a titre d'exemple une 
transformee en cosinus discrete a 1* image d' entree et 

5 un circuit PARE determine les parametres Pe de 1 • image 
d* entree le qui seront utiles a 1' Evaluation a realiser. 
De meme, un circuit TCDS applique a titre d'exemple une 
transformee en cosinus discrete a 1 ' image de sortie I3 et 
un circuit PARS determine les 'parametres Pg de 1 ' image de 

10 sortie Ig qui seront utiles a 1' evaluation a realiser. 

Un circuit de comparaison COMP effectue des 
comparaisons entre les parametres d' evaluation Pe et Ps 
et produit en sortie des resultats d' evaluation R. Les 
parametres Pe sont transmis vers la sortie par le systeme 

15 TRANS ou independamment de celui-ci. Alternativement , il 
est possible de faire "remonter" les parametres Ps vers 
1* entree afin qu'ils y soient compares avec les 
parametres Pe, 

On remarquera egalement qu'il est possible 

20 d'utiliser les coefficients TCD disponibles au niveau du 
codeur COD et/ou du decodeur DEC, mais alors 1' evaluation 
ne tiendra pas compte des defauts du codeur COD et/ou du 
decodeur DEC. 

La transformee en cosinus discrete bi- 

25 directionnelle TCD fait correspondre a un bloc de N x N 
pixels (en general .V = 8) de 1 ' espace pixel un bloc de N 
X N coefficients du domaine transforme. Ces derniers 
expriment la representation frequent ielle de 1* amplitude 
de la luminance des pixels - 

30 La transformee TCD presente des coefficients 

r (u,v) dans lesquels : 

u^v = index de la frequence spatiale 
norizontale et verticale^ 

x,y = index de position horizontale et 

35 vertlcale du pixel dans le bloc. 
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- f{x,y) represente la I'uminance du pixel au 
point (x,y) dans 1 'image, 

- F(u,v) represente le coefficient tranforme 
au point (u,v) du domaine transforme, 

5 - avec u,v,x,y variant de O a N-1 (N = 8). 

Les parametres exploites pour representor les 
degradations contenues dans les images numeriques codees 
notamment en MPEG2 utilisent cette transformation. 

Une information sp'atiale est fournie a partir 
10 du domaine transforme par le caicul de I'activite 
spatiale globale SA de 1' image. Ce parametre informe 
precis^ment sur les details dans le signal image (I,) 
indexe a 1' instant t. (I^) contient (HxW) blocs de pixels. 
La taille du bloc f„„,(x,y) est 8x8 pixels (voir fig.l). 

Dans un premier temps la trans formee TCD par 
blocs est appliquee a 1 • image entierement, ou 
partiellement (sur une zone d'interet ou fen^tre 
d- analyse). Le premier facteur representatif du contenu 
de 1' image, I'activite spatiale, est obtenu k partir de 
I'image tranformee TCD (I,). Un bloc image tres actif 
spatialement a une representation consequente dans le 
domaine transforme. Ses coefficients ont une amplitude et 
une distribution spectrale importantes. 

Deux parametres sont introduits, I'activite 
spatiale de bloc (bsa) et I'activite spatiale d'image 
SA(IJ. lis representent respecrivement I'activite 
spatiale intra-bloc et intra-image. 

Le contenu analyse est celui de chaque bloc de 
pixels f™,(x,y)de position a n,m dans i'image (I.). 

L'indicateur (bsa^,„) pour le bloc n,m tient 
compte des composantes alternatives du bloc de 
coefficients transformes F„,„(u,v) . 

on a : = M/-.«(--->-)) = [i;i:(F,.„C/.y))^l -{/=;^(0.0)y 

•^5 Vi=o j=o J 
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SA(It) est 1' expression du contenu global de 



1* image (ou de la zone d' analyse). 
iTioyenne des acLivites spatiales 
concernes . 

1 



//-I 



:i est calcule par une 
(bsap.^ja) des blocs 
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OU (H X W) est le nombre de blocs de pixels de 
1' image ou de la fenetre d' analyse. 

L'indicateur d'activite spatiale global SA{1^) 
permet de qualifier la richesse de 1' image. Ce parametre 
est affecte par le codage a reduction de debit er par les 
erreurs de transmission. La comparaison du meme parametre 
calcule a 1' entree et a la sortie d'un systeme video 
numerique contribue a 1 ' identification des degradations 
generees par ce systeme. Ce parametre est representatif 
des degradations generees par le systeme, mais sans 
1' inconvenient d'une compensation entre les defauts. En 
effet, d'une part le calcul realise par blocs elimine les 
effets de bloc er d' autre part les pertes d' information 
et les faux contours se traduisent dans la tranformee par 
des effets de parte des hautes frequences, et ces deux 
phenomenes agissenr done dans ie meme sens. 

L' evolution temporelle de SA(I^) represente 
les evenements apparus dans la sequence video. Par 
ailleurs, ce parametre est affecte par le codage a 
reduction de debit et par les erreurs de transmission. La 
comparaison du meme parametre calcule a 1" entree et a la 
sortie d'un systeme video numerique est effectuee sur une 
image I^. indexee a 1* instant t. Elle contribue a 
1' identification des degradations generees par ce 
systeme. La formulation mathematique utilisee pour la 
comoaraison est : 

log_rarzo(/J= logio 
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Le parametre definitif pour, la mesure de la 
degradation liee a la parte d ' informations est ie 
parametre MLR : 

MLR^, = ^^aximui%[log_ratio ( I^) ] 

Get indicateur est calcule sur une sequence de 
M images. Cette option permet d' avoir une information 
reduite concernant une periode de temps. Elle permet 
d* avoir une note de qualite sur : T = M x 1/25 seconde, 
ou plus generalement T = M/NI, NX designant le nombre 
d" images par seconde. 

La procedure de calcul de la perte 
d' information se presenre par exemple comme suit : 

a) selection d'une image a 1' entree du 

systeme, 

b) calcul de SA{1^)^ pour cette image, 

c) indexation de cette image avec sa date et 
son parametre associe S^^, 

d) transmission de parametres de reference 
vers la sortie du systeme, 

f) calcul de SA(l^)^ pour cette image, 

g) calcul du log.ratio(I^) pour chaque image 

T 

h) calcul de la no-ce de qualite MLR pour 
chaque sequence de M images. 

On notera que le proced6 pourrait egalement 
etre mis en oeuvre en faisant "remonter" vers 1' entree 
les valeurs de SA calculees pour chaque image de sortie. 

Un avantage du precede est que chaque image a 
1' entree et a la sortie est caracterise par la valeur du 
parametre SA, ce qui signifie qu'il n'est besoin de 
transmettre en tant que parametre Pe qu'une seuie donnee 
par image pour caracteriser les degradations (autres que 
I'effet de bloc) par opposition a des precedes de I'art 
anterieur qui necessitent de transmettre un grand nombre 
de donn6es permettant d'effectuer une comparaison par 
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difference pixel a pixel encre les images reference et 
degradee. 

La comparaison effecruee a partir du facteur 
MLR rraduir une mesure de la qualite des images finales 

5 et par consequent: informe sur les performances du systeme 
de codage mis en oeuvre lors de la transmission. 

Afin de valider la pertinence du facteur MLR^ 
nous avons confronte les resultats obtenus aux notes 
obtenues par des panels de personnes representatives lors 

10 de tests subjectifs, Le fort coefficient de correlation 
obtenu entre 1* evaluation objective et subjective nous 
demontre I'efficacite d'un tel parametre. 

Dans les sequences video codees a faible 
debit, la distorsion par effet de bloc se manifeste par 

15 des discontinuites du contenu de 1* image aux limites des 
blocs adjacents. Cet artefact est visible sur les zones 
de faible ou de grande activite spatiale et aussi sur les 
zones de forte ou de faible luminosite. II est du a un 
codage inter-blocs independant et a une quantification 

20 grossiere du contenu de chaque bloc. Le coefficient 
continu F(0,0) de la transformee TCD represente la 
moyenne des valeurs de luminance d'un bloc- Le codage des 
coefficients continus de deux blocs, adjacents 
spatialementr a des valeurs sensiblement differentes est 

25 une des causes de cette degradation, Les autres 
coefficients transformes, qui representent les 
composantes alternatives, traduisent les differentes 
evolutions contenues dans le bloc. Le codage et la 
quantification grossiere de ces coefficients genere 

30 egalement 1* effet de bloc dans les zones des contours - 
Des erreurs dans le flux numerique influent sur la 
compensation de mouvement. Ceci se traduit par un effet 
au niveau des blocs ou des macro-blocs. 

Selon 1' invention, un parametre de mesure de 

35 1* effet de bloc met en oeuvre une comparaison de deux 
parametres calcules de chaque cote du systeme video. 
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L'indicateur propose permet d'6valuer 
1' importance de 1 ' information visuelle erronee apport^e 
par la grille structuree de blocs 8x8 visible dans les 
images en sortie de chaine. 
^ Les bordures des blocs 8x8 ou les 

discontinuites de contours sont considerees comme une 
information perceptible et additionnelle dans 1' image. 
Sachant que le codage est realise de maniere independante 
pour les blocs adjacents, -les pertes d' information 
10 internes aux blocs generent des discontinuites de 
contenu, voire meme 1' apparition de nouveaux contours sur 
les frontieres. 

Un procede adapte a cette degradation selon 
1' invention repose sur le calcul de deux parametres : 
I'activite spatiale globale SA (voir supra) et I'activite 
spatiale globale decalee SA^ (voir fig. 2). 

Le parametre SA reste non affecte par les 
effets de bloc, puisqu'il est realise a partir de ces 
memes blocs et de maniere alignee sur leur grille. Afin 
de souligner la contribution des contours et des defauts 
autour des blocs, on determine I'activite spatiale 
decalee SA^,. Elle permet de mettre en exergue les 
artefacts aux frontieres. A partir de la meme definition 
de SA, SAfl est calculee avec de preference un decalage de 
N/j = 4 pixels horizontaleraent et de N/2 = 4 pixels 
verticalement. On remarquera que I'activite spatiale 
decalee SA^ peut etre calculee en effectuant un decalage 
de 1 a N-1 pixels (et de preference N/, pixels) seulement 
dans la direction horizontale, ou seulement ' dans la 
direction verticale, de maniere a evaluer I'effet de bloc 
horizontal ou bien I'effet de bloc vertical. 

Les deux parametres sont issus d'une somme sur 
la zone d'interet ou sur 1 ' image entiere. Une comparaison 
entre les deux indicateurs met en valeur I'apport de 
1' information des contours au bord des blocs. Un nouveau 
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parametre, est done introduit, Sa^^^i^ dans 1' equation 
suivante. 

SA^ (t) 
SA{t) 

^^ratio n'est pas seulemenr sensible aux contours ou 
discontinuites aux bords des blocs mais il I'est 
egaiement au contenu d'origine des blocs. II est par 
consequent necessaire, pour mesurer specif iqueraent la 
degradation introduite par. un systeme video, d'en 
extraire 1 • information apportee par le contenu original 
de 1 • image . 

L'indicateur definitif pour la mesure et la 
detection de I'effet de bloc introduit par un systeme 
video est calcule a partir des valeurs SA^^tio obtenues a 
1' entree et a la sortie du-dit systeme. 

Ces valeurs sont considerees respectivement 
comme valeur du signal reference a 1 ' entree (SAratio^) et 
valeur du signal degrade a la sortie (SAratiOs) • 
L'indicateur de I'effet de bloc BM est defini par 
20 1' equation suivante. 

SAratio s 

Le precede de detection de i'effet de bloc 
25 peut etre mis en oeuvre par exemple comme suit : 

a) selection d'une image a 1' entree du 

systeme, 

b) calcul de SA(I.Je. SA^(It)e et SA,a^^o(It)e 
pour cette image, 

indexation de cette image avec sa date et 
son parametre associe, 

d) transmission de donnees de reference vers 
la sortie du systeme, 

e) identification de la meme image indexee a 
35 la sortie du systeme^ 
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f) calcul de SAdJs, SA^(1.)^ et SA.3^io{It)s 
pour la meme image a la sorrie du sysreme, 

g) calcul du BMs^ pour chaque image I^^ 

h) calcul de la moyenne de BMg^ pour chaque 
5 sequence de W images (avec par exemple M' = M) . Cette 

note qualifie la sequence pour la degradation en 
question. 

La figure 3 demonrre 1 ' importance decroissante 
de BMsA pour 4, 6 et respect ivement 9Mbits/s, Cette 

10 remarque est valable specialement dans la derniere partie 
de la sequence video. En effet, pendant cette periode le 
contenu spatial et temporel de la sequence choisie evolue 
brutalement. Ceci provoque une baisse de la qualite des 
images finales due a une reduction de debit. Le codage 

15 par quantification des coefficients TCD provoque 
1' apparition de I'effet de bloc, et ce d'autant plus que 
le codage est grossier. 

Les degradations dans les sequences d' images 
numeriques codees ont des formes spatiales ou 

20 temporelles. En effet, des erreurs de codage ou de 
transmission peuvent affecter le contenu spatial : 
contours, objets et structures contenus dans le signal. 
Cependan-, elles peuvern: egalement entacher 1' evolution 
temporelle de la sequence video : saccades, gel d' images 

25 ou deplacements errones d' objets. 

Afin d' identifier les degradations temporelles 
qui alterent les sequences video, nous avons defini un 
parametre supplementaire appele "Activite Temporelle" TA. 
II traduir les evolutions du contenu des sequences 

30 d' image dans le temps. II est sensible aux differents 
evenements qui peuvent survenir et affecter le signal. 

L'indicateur spatial SA(I^) traduit le contenu 
de 1' image courante dans la sequence a un instant t. Les 
variations temporelles de SA(I.) representent les 

35 evolutions temporelles du contenu de 1' image et done 
elles traduisent les changements survenus (voir fig. 4a). 
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Le signal mono-dimensionnel SA(I^) d'activice spatiale 
globaie de 1 * image esr proporrionnel a la richesse du 
contenu. Dans le cas mono-dimensionnel, la TCD est une 
transformation qui fait correspondre a un ensemble de M 
points /(x) une serie de M frequences. Elle se d6finit 
par 

1 1 



c(u) = xNl si u = 0 . 
c(u) =1 si u 96 0 

F(0)= J^c(0)X/(-v) 



Sur route fenetre de M echanrillons temporels 
15 du signal SA(I^) est appliquee une transformee en cosinus 
discrete TCD. Cette nouvelle representation traduit les 
frequences temporelles du signal et par consequent 
I'activite de la sequence d' image. La premiere raie F(0) 
est la composante continue, les autres raies F(u) avec 
20 u^o sont les composantes alternatives. 

La composante continue du spectre renseigne 
sur la valeur moyenne de la richesse du contenu sur M 
images. Les frequences alternatives informenr sur les 
variations temporelles du contenu. Seule 1* amplitude des 
25 raies est prise en consideration. 

TCD- = TCD {SA{I,)) = I^Fq{u) u = 0.. M-1 

Le spectre TCD^ (voir fig. 4c) du signal SA 
30 (voir fig. 4b) informe sur les evolutions temporelles de 
la sequence video et par consequent traduit toute 
degradation modifiant ces evolutions. F(0) par definition 
represente la moyenne du signal SA(I^) pendant les M 
images de sortie considerees dans la sequence. Les 
35 amplitudes des composantes alternatives F(u)u^o traduisent 
la maniere avec laquelle change le contenu de la video : 
variations rapides (clips video) , modifications lentes 
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("travelling"), plans fixes, etc... Afin d'exploiter les 
oroprietes de 1 'analyse dans le domaine transforme, un 
parametre TA a ete defini. II est calcule pour un groupe 
de M images . ^ 

TA indique 1' amplitude des fluctuations de 
I'activite SA dans la sequence. II traduit tout evenement 
temporal de la sequence et pa-r consequent : les coupures 
et les gels d' images. TA est compart a un seuil donne Sj 
etabli empiriquement , en dessous duquel on considere 
qu'on se trouve en presence d'une coupure ou d'un gel 
d ' images . 

Le procede de detection des coupures et des 
gels d' image peut etre mis en oeuvre par exemple comme 

suit : 

a) selection d'un groupe d' images, 

b) calcul de SA(IJ pour chaque image, 

c) calcul de la TCD^ et de TA, 

d) - si TA<Si et si FCOXSj, alors le groupe 
d' images est une succession de la meme image uniforme (en 
pratique coupure noire) . 

- si TA<Si et si F(0)>S2, alors le groupe 
d' images est une succession de la derniere image non 
erronee gelee par le decodeur. 

- si TA>Si, retour a a. 

Le phenomene de gel d' images est une 
degradation caracterist ique qui apparalt lors d'erreurs 
sur la chaine complete de distribution du signal de 
television. Le parametre d'activite temporelle TA est 
applicable dans ce cas de figure. En effet, il detecte 
cet evenement qui est lie A la strategic du decodeur. Le 
gel d' images prend deux formes : 

1- restitution des dernieres donnees images 
non corrompues. 
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2. visualisation d'une image uniforme (souvent 
noire) , et ceci tant que les donnees . regues sont 
inexploitables . 

La famille de codage MPEG introduit la notion 
5 de groupe d' images ("Group of Pictures GOP"). Cette 
structure de la sequence est liee a une syntaxe 
definissant les algorithmes predictifs a reduction de 
debit qui sont appliques a chaque image (images codees 
sans prediction : Intra - I, images predictives 

10 unidirectionnelles P, images predictives Bidirectionelles 
B) . La particularity de chaque image influe sur 1" aspect 
de 1' impact des differentes erreurs qui I'affectent. 

L*indicateur defini ici traduit les evenements 
temporels qui surviennent durant la sequence video. II 

15 permet par consequent d' identifier et de comptabiliser 
les ruptures temporelles des sequences. Par la meme 
occasion il represente un parametre de qualification de 
la qualite du service. 

L' invention ne se limite pas aux modes de 

20 realisation decrits. En particulier, I'effet de bloc 
pourrait etre evalue de naniere intrinseque sur 1 ' image 
de sortie. Le precede est alors caracterise en ce qu'il 
comporre les etapes suivantes : 

di) determiner l*acrivite spatiale globale SA^,2 

25 de 1* image de sortie dans une fenetre d' analyse decalee 
presentant des blocs de points image ou pixels qui sont 
decales d'au moins un pixel dans la direction des lignes 
de 1' image et/ou dans une direction perpendicuiaire a la 
direction des lignes de 1 ' image par rapport auxdits blocs 

30 de pixels de ladi^e fenetre d' analyse, cette 
determination mettant en oeuvre les sous-etapes suivantes: 
i) appliquer ladite transformee par blocs sur 
chaque bloc decale {n,m) dudit ensemble de blocs de 
pixels de la fenetre d' analyse decalee pour determiner 

35 les coefficients transformes F^jn^m (i,j). 
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ii) determiner a partir des coefficients 
transformes F^n^m (i.j) des blocs decales, I'activite 
spatiale bsa^^ de chaque bloc dudit ensemble de blocs 
decales, 

5 iii) determiner I'activite spatiale globale 

SAd2 de 1' ensemble des blocs decales constituant la 
fenetre d' analyse decalee^ 

ei) comparer I'activite spatiale globale SA^i2 
de 1' image de sortie dans la fenetre d' analyse decalee et 

10 I'activite spatiale globale SAj de 1' image de sortie dans 
la fenetre d' analyse. 

L' invention peut etre mise en oeuvre a I'aide 
de tranformees par blocs lineaires et orthogonales telles 
que la tranformee en cosinus discrete TCD, la trans formee 

15 en sinus discrete, la transformee de Fourier, la 
tranformee de Haar, la transformee d'Hadamar, la 
transformee de Slant ou bien encore une transformee en 
ondelettes (DWT) dont les coefficients ont ete 
reorganises pour obtenir les blocs transformes. 
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REVENDICATXONS 

1. Procede d' evaluation de la degradation 
d'une image video codee par blocs de points images ou 
pixels, cette degradation etant generee par un systeme de 
5 codage et/ou de stockage et/ou de transmission, 
produisant une image de sortie a partir d*une image 
d' entree, caracterise en ce qu'il comporte les etapes 
suivantes : 

a) selectionner une image d' entree et 
10 determiner une activite spatiale (SA) de 1 • image d* entree 

dans une fenetre d' analyse representant au moins une 
partie de 1 ' image qui presente un ensemble de dits blocs 
de points images ou pixels, cette determination mettant 
en oeuvre les sous-etapes suivantes : 
15 i) determiner pour chaque bloc (n,m) du dit 

ensemble de blocs de pixels, les coefficients transformes 
Fn,m {i,y) par une transformee par blocs selon les 
dits blocs 

ii ) determiner a partir des coefficients 
20 transformes Fn,m (i,y) 1* activite spatiale bsa de chaque 

bloc dudit ensemble de blocs de pixels 

iii ) determiner a parrir de 1* activite 
spatiale de bsa de chaque bloc 1' activite spatiale 
globale (SA^) de 1' ensemble des blocs de pixels 

25 constituant la fenetre d' analyse 

b) selectionner 1' image de sortie 
correspondant a 1 ' image d" entree et determiner ladite 
activite spatiale globale (SA2) de 1 ' image de sortie dans 
ladite fenetre d' analyse, en mettant en oeuvre pour 

30 1' image de sortie les sous-etapes a i) a a iii) 

c) comparer 1' activite spatiale globale (SAj) 
de 1' image de sortie dans la fenetre d* analyse et 
1* activite spatiale globale (SAO de 1 ' image d' entree dans 
la fenetre d' analyse. 
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2. Procede selon la revendicarion 1 
caracrerise en ce que I'activite spatiale bsa d'un bloc 
(n,m) est determinee par ia formule suivante : 

V«=o/=o J 

ou k(i,j) est un coefficient constant de la ponderation 
pour la composante i,j utilises et p une constante. 

3. Procede selon la revendication 2 
10 caracterise en ce que I'activite spatiale SA de 

1' ensemble des blocs peut #tre alors determinee par la 
formule suivante : 

(HxW) representant le nombre de blocs de 
pixels dans la fenetre d' analyse. 

4- Procede selon la revendication 3 
caracterise en ce que ladite comparaison (etape c) est 
effectuee a I'aide d'un parametre LR defini de la maniere 
suivante : 

I-R = gj Cfi(SAi/SA2) ] avec i,j6 {1,2} 
avec fi(x,y)=x-y ou f2(x,y)=x/y 
et gi(x)=100lx| ou g2(x)=iOO |log(|x|)( 

5. Procede selon la revendication 4 
caracterise en ce qu'ii „,et en ceuvre le calcul du 
parametre LR sur une sequence de M images d' entree et de 
M xmages de sorties correspondants, et en ce qu'il met en 
oeuvre le calcul d'un parametre d' evaluation ou note de 
qualite MLR defini de la maniere suivante : 

MLR = Maximumv (LR) 

6. Procede selon une des revendications 
pr^cfedentes caracterise en ce qu'il comnorte la 
determination d'un indicateur d'activite temporelle TA 
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pour un groupe de M images, determine de la maniere 



suivante 



5 Fs;^(u), pour u variant de 0 a designant 

les M coefficients d'une transform6e par blocs appliquee 
a M echantillons temporels de I'activite spatiale SA. 

7. Proced6 selon la revendication 6 
caracterise en ce qu'il compofte une etape de comparaison 

10 de 1 ' indicateur d'activite temporelle TA a un seuil donne 
Si, de maniere a determiner la presence d*une coupure ou 
d'un gel d* image. 

8. Procede selon une des revendications 
precedentes caracterise en ce qu'il comporte, en vue de 

15 la determination de maniere intrinseque de l*effet de 
bloc/ les etapes suivantes : 

di) determiner I'activite spatiale globale SA^2 
de 1' image de sortie dans une fenetre d* analyse decalee 
presentant des blocs de points image ou pixels qui sont 

20 decales d'au moins un pixel dans la direction des lignes 
de 1 • image et/ou dans une direction perpendiculaire a la 
direction des lignes de 1 ' image par rapport auxdits blocs 
de pixels de ladite fenetre d' analyse, cette 
determination mettant en oeuvre les sous-etapes suivantes 

25 i) appliquer ladite transformee par blocs sur 

chaque bloc decale (n,m) dudit ensemble de blocs de 
pixels de la fenetre d' analyse decalee pour determiner 
les coeeficients transformes F^n^m 

ii) determiner a partir des coefficients 
30 transformes F^n^m (i,j) des blocs decales, I'activite 

spatiale bsa^, de chaque bloc dudit ensemble de blocs 
decales, 

iii) determiner I'activite spatiale globale 
SAd2 de 1' ensemble des blocs decales constituant la 

35 fenetre d' analyse decalee, 
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ei) comparer I'acrivite spatiale globale SA^j 
de 1' image ce sortie dans la fenetre d' analyse decalee et 
I'activite spatiale globale SA, de 1 ' image de sorrie dans 
la fenetre d' analyse. 
5 9. Precede selon une des revendications 

precedentes caracterise en ce qu'il comporte^ en vue de 
la determination de I'effet de bloc, les etapes. suivantes 

d2) determiner I'activite spatiale globale SA^^i 
de 1 • image d* entree dans une" fenetre d' analyse decalee 

10 presentant des blocs de points image a pixels qui sont 
decales d'au moins un pixel dans la direction des lignes 
de 1' image et/ou dans une direction perpendiculaire a la 
direction des lignes de 1' image par rapport auxdits blocs 
de pixels de ladite fenetre d* analyse |. cette 

15 determination mettant en oeuvre les sous-etapes suivantes 

i) appliquer une trans formee par blocs sur 
chaque bloc decale (n,m) dudit ensemble de blocs de 
pixels de la fenetre d* analyse d6calee pour determiner 
les coeeficients transformes F^n^m (i,y)^ 

20 ii) determiner a partir des coefficients 

transformes F^n^m (i,y) des blocs decales, I'activite 
spatiale bsa^j de chaque bloc dudit ensemble de blocs 
decales^ 

iii) determiner I'activire spatiale globale 
25 3A^^ de 1* ensemble des blocs decales constituant la 
fenetre d» analyse decalee, 

e) determiner I'activite spatiale Sh^2 de 
1* image de sortie dans ladite fenetre d' analyse decalee 
en mettant en oeuvre les sous-etapes d^ i) ad iii) 
30 f) comparer d'une part I'activite spatiale 

SA^2 cie 1' image de sortie dans la fenetre d* analyse 
decalee et I'activite spatiale SAj de 1 ' image de sortie 
dans la fenetre d' analyse, pour evaluer I'effet de bloc 
dans 1' image de sortie et d' autre part I'activite 
35 spatiale SA^^ de 1 ' image d'eniiree dans la fenetre 
d' analyse decalee et I'activite spatiale SA^ de 1 ' image 
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d' entree dans la fenetre d' analyse pour evaluer I'effet 
de bloc dans 1' image d' entree. 

10. Procede selon la revendication 9 
caracterise en ce que l'6tape f met en oeuvre la 

5 determination d'un indicateur BM d'effet de blocr selon 
la formule suivante : 

BM = X 100 

10 

11. Procede selon une des revendications 
precedentes caracterise en ce que ladite transformee par 
blocs est une transformee lin6aire et orthogonale^ 
notamment une transformee en cosinus discrete TCD ou une 

15 transformee en sinus discrete ou une transformee de 
Fourier ou de Haar ou de Hadamar ou de Slant. 

12. Procede selon une des revendications 1 a 
10 caracterise en ce que ladite transformee par blocs est 
une transformee en ondelettes (DWT) dont les coefficients 

20 ont ete reorganises pour obtenir des blocs transformes. 
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